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Abstract 

A model of two interacting one-dimensional fermion systems ("Luttinger 

liquids") coupled by single-particle hopping is investigated. Bosonization 

allows a number of exact statements to be made. In particular, for forward 

scattering only, the model contains two massless boson sectors and an Ising 

type critical sector. For general interactions, there is a spin excitation gap 

and either s- or d-type pairing fluctuations dominate. It is shown that the 

same behavior is also found for strong interactions. A possible scenario for 

the crossover to a Fermi liquid in a many chain system is discussed. 
74.20.Hi, 75.10.Lp, 71.45.Lr 
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The properties of a strictly one-dimensional interacting fermion system are by now rather 
well understood. |I],|2| The typical phenomenology (called "Luttinger liquid" 0) is charac- 
terized by a separation of the dynamics of spin and charge and by interaction-dependent 
power laws in many correlation functions, and is thus quite different from Fermi liquid be- 
havior familiar from higher-dimensional systems. On the other hand, the effects of coupling 
between parallel chains, present in any real gwasi-one-dimensional system, are still a subject 
of debate. Considerable effort has been devoted to the understanding of the properties 
of many coupled chains, || however, it is in many respects unclear how to connect these 
results to the strictly one-dimensional case. A possible bridge between the single and many 
chain cases are two (and possibly three, four, etc.) coupled chains. The two-chain case 
is also of relevance for experiments on Sr 2 Cu 4 6 , (VO) 2 P20 7 , [§] and possibly the blue 
bronzes U (in this last case three-dimensional phonons certainly play an important role). 

The two-chain model has thus attracted considerable interest, both analytically [pT0| p~3f| 
and numerically. |I4|-|T6[| Nevertheless, there is little general information on the low-lying 
excitations or on the possible ground state phases. In the present paper I investigate this 
problem for small interchain hopping integral and small intrachain interaction, but with their 
relative size left arbitrary. Using the standard bosonization procedure, a rather complete 
picture of the different possible phases and the excitation spectrum will emerge. It will 
further be shown that the low-energy properties found for weak interactions also exist in 
the strong-interaction limit, suggesting that weak and strong interaction are in the same 
phase of the coupled chain model. 

The model I consider is given by the Hamiltonian 



a Fermi velocity vp and forward and backward scattering interaction g 2 and gi, t± is the 
interchain hopping amplitude, and ip rs i is the fermion field operator for right (r = +) or left 
(r = — ) going particles of spin s on chain i. To start, I neglect the backward scattering 
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(1) 



Here H 12 are the (identical) Hamiltonians of the two chains, each characterized by 



g%. The following analysis is then initially identical to that of ref. | 10| . The Hamiltonian 
is transformed by the following steps: (i) introduce bonding and antibonding operators via 
iprso = (VVsi + VVs2)/v2, tprs-n- = (VVsi — VVs2) /\/2; (ii) introduce charge and spin boson fields 
<Pp,cr;0,TT corresponding to the 0- and tt- fermions, following the standard procedure; (iii) form 
the linear combinations (p u ± = ((p u o ± l/7r )/v / 2 {y = p,cr). The noninteracting Hamiltonian 
(including t±) then takes the form 



H o = Yl dx 



2 

v=p,a 



7T 



(2) 



where U ua is the momentum density conjugate to (j) ua , and the interaction is 

1 



*mt,2 
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,(2) 
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7 =± 



+ / dx cos 26 P - ( cos 2 </><r- + cos 20*- ) • (3) 

2(7T«J / 7 

Here a is a short distance cutoff, d x 8^ = tiTL^,, g\ = <7oo o +7fl , 07nrO> an< ^ ^ use ^ ne notations 
of ref. [pT| : p^L * s ^ ne coupling constant for an interaction scattering two particles from 
states (a, b) into (cf, c). Initially, all the g's in eq. (|3|) equal g%, but renormalization will give 
rise to differences. At energy scales higher then t± an additional process of type g^ty* a l so 
exists and is responsible for the fact that gi is not renormalized in the purely one-dimensional 
problem t± = (this process also only involves the p_ and cr_ fields). At energies below t± 

(2) 

the ^ottOtt - process becomes however forbidden due to energy and momentum conservation, 
and eq.(^) is then indeed the full forward scattering Hamiltonian. 

One now can notice that the p + and a + parts of the Hamiltonian remain bilinear, and 
the corresponding fields are thus massless. On the other hand, there are nontrivial inter- 
action effects for the coupled p_ and <r_ fields: one finds coupled Kosterlitz-Thouless type 
renormalization group equations for g^ n7r and g_ . [[iI|[DJ For the initial conditions appro- 
priate here, these equations always scale to strong coupling, and the standard interpretation 
|I(J then is that there is a gap A « t± exp(— n 2 VF/\g2\) f° r both the p_ and <r_ degrees of 
freedom. 
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That things are actually a bit more subtle can be seen noting that the cr_ part of the 
Hamiltonian is the continuum transfer matrix of a two-dimensional classical XY model with 



twofold anisotropy field cos20 cr _ (the XY spins then are (S = cos</> CT _, sin0 CT _)). [p~8|JT9 
This model has Ising type symmetry, with order parameter sin , and the symmetry of the 
Hamiltonian under the duality transformation <-» Q a _ implies that the model is critical. 
The duality symmetry is related to the fact that the left- and right-going fermions are 
independently invariant under spin rotation, i.e. there is a chiral 577(2) x 577(2) symmetry 
in the fermionic model. 

What are the physical properties of the pure forward scattering model? First, there 
are massless modes in the p + and a + channels, giving a total specific heat C(T) = 
(tvT/3)(1/u p+ + l/u a+ + l/2vp), where the total charge and spin velocities are given by 



u 



2 

P+ 



{92 /tt) 2 and u a+ = Vp, and the factor 1/2 in the last term comes from the 



Ising critical behavior (with central charge c = 1/2 p0[). The compressibility is determined 
by the p + modes only and given by n^ 1 = 7rp 2 ) u p+ /4:K, where po is the equilibrium particle 
density and K 2 = {'KVp — g^)/ (^v F + g 2 ). Similarly, the (Drude) weight of the zero-frequency 
peak in the conductivity is cr = 4u p+ K. As in the one-chain case, |J these relations can in 
particular be used to determine the coefficient K which determines power laws of different 
correlations functions. 

Naturally, the present model does not have broken symmetry ground states, but as in 
the one-chain case there are divergent susceptibilities of different types, indicating incipient 
instabilities. I first consider g 2 > 0. To obtain the long- wavelength (low-energy) asymptotics 
of correlation functions one has to analyze the consequences of the nonlinear term in eq. @ 
which scales to strong coupling (gQ 0nn — > 00). A semiclassical treatment is appropriate, and 
then the energy is minimized by 6 p - = (there are different degenerate solutions which 
all lead to identical physical results). Following standard arguments [|1|] long-range order 
of the 9p- field implies exponentially decaying p _ correlations. On the other hand, from 
the Ising analogy for the <r_ sector correlations of the order parameter sin and its dual 
sin 6 a _ then decay as r _1//4 whereas correlation of the non-ordering cos (f) a _ and cos 9 a ^ fields 



decay exponentially. These points have not been appreciated in previous work on this model. 
Consider now for example charge density oscillations which are out of phase between the two 
chains, described by the operator Ocdw-k ~ e l ^ p+ sin0o.+ sin0 CT _. From the massless modes 
the CDWjr correlations then decay as r -( 3 + 2/ 0/ 4 ; giving rise to a susceptibility diverging 
as T( 2X - 5 )/ 4 . The analogous spin (SDW^) correlations obey the same power law, whereas 
in-phase correlations decay exponentially. 

Similar considerations apply to BCS type instabilities. It turns out that long-range 
correlations exists for the pairing operator 

S Cd = ^ S (^- - Sl 0^+, S ,0 - 4>--s,k4>+,s,it) (4) 
s 

and its triplet analogue. It seems appropriate to call this form "d-wave" because pairing 
amplitudes of the "transverse" modes and ir intervene with opposite sign. The bosonic 
form of this operator is given by the same form as Ocdw-k, with </> p+ — > 9 p+ , 6> CT _ — > The 
corresponding susceptibilities diverge like T^ 2// ^ -5 ^ 4 . Because for g 2 > one has K < 1 this 
divergence is subdominant compared to the CDW^ and SDW,,- ones. It may seem surprising 
that the exponents do not tend to zero as g 2 — > 0, however one should notice that the power 
laws are valid for T < Aq, and because A — + for g 2 — > there is a nontrivial crossover 
in the noninteracting limit. In all other pairing correlations, "s-wave" superconductivity 
in particular (a plus instead of the minus sign in eq.(f|), the leading divergent terms cancel 
and one therefore has exponential decay of correlation functions and finite susceptibilities 
as T — > 0. 

(2) 

For negative g 2 the picture changes quite drastically, because now scaling goes to gQ 07r7T — > 
— oo, and consequently the Ising order parameter is cos</> CT _. Now K > 1, and the dominant 
divergent susceptibility is then easily found to be standard s-wave superconductivity, with 
exponent {2/K — 5)/4. The subdominant divergence occurs for orbital antiferromagnetic 
operators [^TJ of the form ip\ S7 ,ip_ s() — ipl_ s0 ip_ S7r and its triplet analogue (the spin nematic). 

Consider now the backscattering interaction g\. I will only treat the repulsive case 
gi > 0. In a purely one-dimensional system this then scales to zero as gi(E) = g\/(l + 
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gi/ (tcvf) \tl(vf/ Ea)) when the running cutoff E goes to zero. In the coupled chain problem, 
the one-dimensional scaling breaks down for E ~ t±. For small t± the effective g\ = gi(t±) 
will then indeed be a perturbation. Simultaneously, g 2 is renormalized to g 2 = <?2 — gi/2 + 
g{/2. The backscattering Hamiltonian takes the form 
2g* f 

Hmi = 7 / rfx{cos20 cr+ (cos 2^ p _ + cos20 cr _ + cos26> CT _) — cos2#„_ cos2# CT _} 

(27ra)^ J 

(5) 



•— / dx 



^{{dAp+f + (d x <p a+ ) 2 } - U 2 p+ - n 2 + 

First, the 6 p -—6 a - interaction now breaks the self-duality of the fields. As cos2# p _ 
already has a nonzero expectation value from the g 2 interaction, one now also finds a gap in 
the cr_ modes, of order A a = (gl/g^Ao. In the Ising model language, this corresponds to a 
deviation from criticality, long-range order and exponentially decaying sin correlations. 
Secondly, the forward scattering interaction also leads to a nonzero expectation value of 
cos 26 p - + cos 20 CT _ + cos 26 l -_ which by spin rotation invariance has to be positive. The 
leading order effect of the first term in eq.(|3p then is to open a gap also in the a + degrees 
of freedom, given, up to numerical factors, by A CT . In the presence of the backscattering 
interaction there thus is a gap in all the magnetic excitations. 

In correlation functions, to leading order one now replaces <p a+ by its classical value tt/2. 
One then finds for #2 > a decay of the SCd correlations as r~ 1 ^ 2K , giving rise to a divergence 
of the corresponding susceptibility as y 1//2K ~ 2 , where now K 2 = (7cvF—g2+gi/2)/(nv F +g 2 — 
<7i/2). On the other hand the CDW Z and SDW„- operators contains the Ising disorder field, 
and therefore these correlations decay exponentially. A divergent density response exists for 
correlations of the form {[Ocdw^)] 2 [Ocdwtt{^)} 2 ) oc cos(2(kp + kF n )r)r~ 2K , |L3| because 



here the operator sin 2 ^- « 1/2 appears. Perturbative and symmetry arguments show 
that the same contribution also exists in the density correlations: (n(r)n(O)) oc cos(2(A;po + 
kp n )r)r~ 2K , in analogy with the Akp oscillations of a single chain. However at least for weak 
interactions (K — > 1) the corresponding susceptibility is much weaker than the SCd pairing, 
i.e. the two-chain model has predominant pairing fluctuations even for purely repulsive 



interactions. |22| In the regime of negative gi the leading divergent susceptibility is s-wave 



superconductivity with exponent 1/2K — 2. The precise boundary between the two regimes 
can be determined from the scaling equations of ref. [XT] and is given by g\ = 2g 2 . The 
triplet susceptibilities (spin density wave or triplet superconductivity) are suppressed by the 
spin gap. The spin gap gives rise to "anomalous flux quantization" , and there is also a 
gap for single-particle excitations. 

The power laws discussed above apply in the temperature region below Ao- In the 
intermediate region A < T < t± the g^ n7T term in eq.(f|) has little effect, and one then can 
obtain the temperature dependence of different correlation functions from a purely bilinear 
Hamiltonian. For example, for CDWo susceptibilities one finds a power law T^ -1 ^ 2 , whereas 
in the one-dimensional region T > t± one has a behavior as T K ~ l . The important point 
here is that in the intermediate region the interaction dependent exponent is smaller than 
in the high-temperature region, i.e. below t± the system behaves more closely like a Fermi 
liquid than at high temperatures. 

Let me now briefly consider the strongly interacting case. For sufficiently strong in- 
trachain interactions, i.e. small parameter K p of the individual chains, single-particle 
hopping is renormalized to zero, however simultaneously particle-hole tunneling processes 
appear. []5]|| Introducing 4> v ± = (<p v i ± J/2 )/\/2, where 4> u i,2 are the boson fields of 
the individual chains, for the purely forward scattering case, these terms take the form 
J cos 20 p _ (cos 20o-_ + cos2#o-_). One again has a duality symmetry, </> CT _ <-> 9 a _, and the 
same types of power-law correlations as in the weak-coupling case appear. Introducing now 
intrachain backscattering, the duality symmetry is broken and, again as in the weak-coupling 
case, only SCd correlations (exponent 1/2K P ) and 4kp charge correlations (exponent 2K P ) 
remain. The types of possibly divergent response functions and the scaling relations between 
different exponents are thus identical for weak and strong interaction. This strongly sug- 
gests that this type of behavior actually holds for arbitrary interaction strength. Note that 
the density correlations decay more slowly than the pairing correlations only for K p < 1/2. 
This typically corresponds to rather strong repulsion: for example, in the one-dimensional 
Hubbard model one reaches K p = 1/2 only for infinite repulsion. |J Another interesting 



strong-coupling model is the "t— J ladder". [TjJ Here in the limit of strong interchain ex- 
change a mapping onto an effective single-chain hard core boson model can be made, leading 



again to the same low-energy properties as in the weak-coupling limit. [E3] Recent numerical 



results (nynj confirm this point. 

The exponents K — 1 and (K — l)/2 valid for the single and double chain problems 
suggest that for N chains coupled by near-neighbor interchain hopping one might have an 
anomalous exponent (K — 1)/N at T < t±. To see how such a behavior can possibly arise, 
in analogy to the two-chain case one can go to momentum space in the transverse direction. 
The noninteracting bosonized Hamiltonian then is 



Ho 



TTVp 



k 



TT 2 
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(6) 



Following standard arguments |25| I now only consider forward scattering interactions which 
for states at the Fermi energy are consistent with both energy and momentum conservation. 
The analogue of the first term in eq.(|3]) then is 



int,2 



— I dx 

2 



\{d x <Ppo) 2 

7T Z 



n 



pO 



(7) 



where <p px is the Fourier transform of <ft p k ± with respect to k±. The important point here 
is that only the mode at x = is affected by the interactions. A standard calculation 
then leads to a decay of CDW correlations as r - 2 -( ftr - 1 )/ Ar ) giving rise to a susceptibility 



behaving as T^ K . Similarly, the single particle Green function decays as 



-1-5 



with 



5 = (K + 1/K — 2) / 4N, leading to a singularity of the momentum distribution function as 
\k — kp\ 5 . p6j In the limit of a large number of coupled chains the anomalous exponents now 
vanish, and in particular one recovers a Fermi liquid like momentum distribution function 
in this description. 

Clearly, a number of interactions has been neglected in this argument. First there are 
Cooper type ((k, —k) — > (k', —k')) and possibly nesting interactions, the prototype of which 
is given by the ^qotttt t erm i n ec L-([l)- By analogy with that case I expect these interactions to 
give rise to a gap of order Aq, and to ordered ground states for N —>■ oo. Thus the power laws 



S 



of the preceding paragraphs are valid in the temperature region Aq < T < t±. Moreover, 
there are interactions that involve at least one state not exactly at the Fermi energy. Though 
these interactions cannot directly affect the low-energy physics, they in general will lead to 
renormalizations of g%. The above arguments remain valid if these renormalizations are 
nonsingular. To which extent this is correct is currently under investigation. 

In conclusion, I've investigated the phase diagram and excitation spectrum of two Lut- 
tinger liquids coupled by single-particle hopping, and proposed a possible extension to many 
coupled chains. The conclusions are valid for small hopping amplitude, but the same types 
of divergent responses (d-type superconductivity and 4k p charge density in the case of 
repulsion) occur for both weak and strong interactions, suggesting that this type of be- 
havior is to be found for rather general interactions. The fact that for strong interactions 
interaction interchain hopping renormalizes to zero only affects properties at inter- 

mediate energy scales (above the spin gap). Contrary to the case of a single chain, the pure 
forward-scattering model is found to be a singular line in the phase diagram, with Ising 
type criticality. 

I am grateful to L. Balents, T. Einarsson, M.P.A. Fisher, T. Giamarchi, R. Noack, 
D. Poilblanc, J. P. Pouget, and H. Tsunetsugu for stimulating discussions. Laboratoire de 
Physique des Solides is a Laboratoire Associe au CNRS. 
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